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Hyperfine-induced spin relaxation of a diffusively moving carrier in low dimensions: 
implications for spin transport in organic semiconductors 

V. V. Mkhitaryan and V. V. Dobrovitski 
Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA 

The hyperfine coupling between the spin of a charge carrier and the nuclear spin bath is a pre¬ 
dominant channel for the carrier spin relaxation in many organic semiconductors. We theoretically 
investigate the hyperfine-induced spin relaxation of a carrier performing a random walk on a d- 
dimensional regular lattice, in a transport regime typical for organic semiconductors. We show that 
in d = 1 and d = 2 the time dependence of the space-integrated spin polarization, P(t), is dominated 
by a superexponential decay, crossing over to a stretched exponential tail at long times. The faster 
decay is attributed to multiple self-intersections (returns) of the random walk trajectories, which 
occur more often in lower dimensions. We also show, analytically and numerically, that the returns 
lead to sensitivity of P(t) to external electric and magnetic fields, and this sensitivity strongly de¬ 
pends on dimensionality of the system (d = 1 vs. d = 3). Furthermore, we investigate in detail 
the coordinate dependence of the time-integrated spin polarization, <r(r), which can be probed in 
the spin transport experiments with spin-polarized electrodes. We demonstrate that, while <r(r) is 
essentially exponential, the effect of multiple self-intersections can be identified in transport mea¬ 
surements from the strong dependence of the spin decay length on the external magnetic and electric 
fields. 

PACS numbers: 72.25.Dc, 75.76.+J, 85.75.-d 


I. INTRODUCTION 

Spin dynamics of the charge carriers in organic semi¬ 
conductors have attracted much attention recently^ - — . 
On one hand, the processes which underlie the use 
of these systems in the organic light-emitting diodes 
and organic photovoltaic solar cells, are explicitly spin- 
dependent, so that understanding of the spin dynamics 
in organic semiconductors is of fundamental interest for 
such applications. On the other hand, the long spin life¬ 
times of the carriers in the organic semiconductors make 
them an interesting candidate for prospective spintronic 
applications, and the detailed understanding of the mech¬ 
anisms of the spin relaxation is required. 

Typically, charge transport in organic semiconductors 
occurs via random (incoherent and inelastic) hopping of 
the polarons carrying positive or negative charge between 
localized molecular sites. During the waiting time be¬ 
tween two consecutive hops the carrier spin interacts via 
hyperfine coupling with the spins of the nuclei (mostly, 
hydrogen), which surround the host site. Thus, the spin 
of the carrier waiting for the next hop undergoes rota¬ 
tion around a random axis by a random angle at each 
site. The principal role of this mechanism was experi¬ 
mentally confirmed* and supported by a number of spin 
transporl£ii2 and spin resonanc e 1 ’ 13 ' 14 measurements, as 
well as theoretical worksi 15-18 Another source of spin 
relaxation is the spin-orbit interaction, which leads to 
the spin-flip scattering in the course of a hop, but this 
interaction in organic semiconductors is weak, so the 
hyperfine-induced spin relaxation is likely to be the main 
source of depolarization, although other theoretical^ 
and experimental^*^ studies favor the spin-orbit mech¬ 
anism, and theoretical efforts are made to explain this 
controversy^ 


Spin dynamics of the polarons in low-dimensional sys¬ 
tems is of particular interest: in the polymer-based de¬ 
vices the carriers mostly move along the ID polymer 
chains, while hopping from one polymer chain to another 
happens mostly at the intersections. One-dimensional or¬ 
ganic polymer wires can be prepared and their properties 
can be studied in detail, see e.g. Refs. I21II22I . Also, en¬ 
gineering low-dimensional systems is a promising way to 
design organic materials with large magnetoresistance^ 
which are of much technological interest. 

This motivates the theoretical study, described below, 
of the hyperfine-induced spin relaxation of a carrier dif¬ 
fusing via random walk in d = 1 and d = 2 dimensions. 
We investigate in detail not only the (space-) integrated 
spin polarization P(t) (which has been addressed in some 
previous worksi^* 24- —), but also the time-integrated po¬ 
larization <r(r) at the given point r in space, which can be 
measured in spin-transport experiments: 2-11 ' 27 We show 
that for low-dimensional transport, these two quantities 
are related to each other in a rather non-trivial way. 
Moreover, we analyze the spin decay length Is for d = 1 
and d = 2, and show that it is very sensitive to both elec¬ 
tric and magnetic fields; this is important both for appli¬ 
cations and for the fundamental studies of the transport 
in organic semiconductors. 

The average polaron hopping rate v, corresponding to 
typical mobilities of 10 -8 to 10 -6 cm 2 V -1 s -1 , is of the 
order of 1-100 GHzi 2£ At the same time, the average 
hyperfine-induced spin precession frequency is of order 
100 MHzj&*!£ Therefore, the hopping is much faster than 
the hyperfine precession, and the carrier performs many 
random-walk steps before its spin polarization averages 
out to zero. Thus, the spin relaxation should be sensi¬ 
tive to the statistics of the underlying random walk. The 
particularly important feature of the random walk is the 
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FIG. 1: (Color online) Lateral (a) and vertical (b) spin valves 
with organic active layer (gray), in which the charge carri¬ 
ers are restricted to move in one-dimensional channels along 
the ^-direction. The separation between the injecting and 
detecting electrodes is L. 


frequency of returns of the hopping carrier to the same 
site, i.e. the frequency of the self-intersections of the ran¬ 
dom walk trajectories. When the returns are absent, in 
the so-called transient diffusion regime^, the random hy- 
perfine field acting on the carrier’s spin has no memory: 
the local hyperfine environments at different sites are un¬ 
correlated, and the carrier hops from one site to another 
never coming back. This corresponds to the motional 
narrowing regime of the spin relaxation^ and leads to 
the exponential decay of the space-integrated spin polar¬ 
ization, P(£) oc exp(— t/ts,ti), where ts,t r is the spin re¬ 
laxation time for the transient-diffusion regim e^- 1 - 32 . The 
opposite scenario with frequent returns, known as the 
persistent diffusion, takes place when the carriers move 
in low-dimensional systems. The frequent returns to the 
same site lead to the random hyperfine field with long 
memory (long-time correlations), which suppresses the 
motional narrowing and leads to faster spin relaxation. 
The role of the memory of the hyperfine field has been 
studied for some scenarios, and has been found to change 
the spin relaxation of the spin polarization in d = 1 case 
to P(£) oc exp [— (t/ts,i) 3 ^ 2 ] at short timesi&2£ - — , with 
the decay time Isa depending on the particular model for 
the hyperfine fields. Also, it has been noticed that the 
returns make P(£) very sensitive to the external magnetic 
fields 

However, the detailed knowledge about the carrier’s 
spin relaxation in the case of low-dimensional transport is 
still largely lacking, and our work aims at filling this gap. 
We investigate analytically and numerically the spin re¬ 
laxation for lateral and vertical spin valves (Fig. [I]), with 
the carriers moving along the ID current-carrying chan¬ 
nels. We consider the limit of small current density, when 
each carrier moves independently of others, thus working 
within the single-particle framework. For the lateral spin 
valve, Fig. Oja), the spin carrier hops along the very long 
linear chain, performing an unbounded random walk. For 
the vertical spin valve, Fig. [IJb), the random walk hap¬ 


pens over finite-size chain with reflecting boundaries. In 
both cases we assume that the carrier is injected with the 
spin state “up” at x = 0, and its spin is probed by the 
detector at x = L. 

We study the space-integrated spin polarization P(t) 
and the time-integrated polarization <r(r) at the given 
point r, which can be measured with the detector lead 
at a given location. We demonstrate that these quanti¬ 
ties exhibit remarkable universal features, and that the 
returns in the course of the carrier diffusion play an im¬ 
portant role in this relation. 

First, for the time decay of P(t), we demonstrate the 
evidence of the universal scaling: for different values of 
the hyperfine coupling strength fohf and hopping rate 
z/, the decay of P(t) follows the same curve which de¬ 
pends only on the normalized dimensionless time r = 
(vt)(bh{/v) 4 ^ 3 . This scaling holds not only for short 
times, where P(t) follows the previously known decay law 
P(£) = exp [—(£/£s) 3 / 2 ], but also at long times, where we 
find a previously unnoticed stretclied-exponential decay 
P(£) ~ exp(—out 3 / 4 ) (for both spin valve geometries). 
This scaling holds for both lateral and vertical spin-valve 
geometries, at finite magnetic and electric field. Also, 
we found that P(£), besides the known sensitivity to the 
magnetic field, is also very sensitive to the electric field. 
We have observed the similarly strong effect of the re¬ 
turns for the spin relaxation in d = 2, leading to the 
logarithmic corrections to the standard exponential de¬ 
cay of P(£), and strong sensitivity to the external fields. 

Second, we studied the time-integrated polarization 
<r(r), which is of much importance for the spin-dependent 
transport measurements. We are not aware of any ana¬ 
lytical theory for this quantity, but our numerical stud¬ 
ies reveal unexpected universality in its behavior. For 
persistent diffusion in low dimensions, in contrast to the 
transient diffusion in 3D, the quantity <r(r) is not di¬ 
rectly related to P(£). Our numerical results show that, 
despite the essentially non-exponential decay of P(£), the 
spin transport decay is exponential, <r(r) oc exp(— r/ls), 
with high accuracy, even in the presence of the external 
magnetic and electric fields. However, the resemblance 
to the usual 3D transient-diffusion result is superficial: 
for both d = 1 and d = 2 cases, the dependence of the 
spin decay length Is on the external fields is very strong, 
in contrast to the standard 3D diffusion. 

Our results suggest that the character of the carrier 
diffusion in an organic semiconductor can be studied in 
spin transport experiments, via the field dependence of 
Is, and vice versa, the spin transport measurements in 
low-dimensional organic semiconductors can be used for 
accurate sensing of electric and magnetic fields, and for 
other similar spintronic applications. 

The rest of the paper is organized as follows. In the 
next Section we discuss the formulation of the problem 
and the methods used for analytical and numerical stud¬ 
ies. In Sections lIIIl and llVI we consider the spin relaxation 
for the lateral and the vertical spin-valve geometries, re¬ 
spectively. Section [V] outlines our results on the spin 
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relaxation in d = 2. Details of the analytical calculations 
are presented in two Appendices. 


II. MODEL FOR THE CARRIER SPIN 
RELAXATION 

We consider a carrier hopping between sites, see 
Fig-HIa), which model organic molecules or conjugated 
segments of polymers. Everywhere below, for both d = 1 
and d = 2, we enumerate sites by the integer variable r, 
so that e.g. for ID chain the physical coordinate of the 
site is x = ar , where a is the distance between the sites. 

When the polaron is localized at the site with the 
radius-vector r its spin interacts with N nuclei I r fc (k = 
1 , ..,7V) surrounding the given site. Below we assume 
that all nuclei have spin 1/2, since the protons in many 
organic semiconductors are the most abundant species 
with the largest nuclear magnetic moment. The Hamilto¬ 
nian governing the spin dynamics of the carrier localized 
at the site r is: 

N 

H T = BS Z + Sj2^rklrk, (1) 

fe= 1 

where a r t is the hyperfine coupling constant between the 
carrier spin and nuclear spin I r fc, and B is the Larmor fre¬ 
quency of the carrier spin S in an external magnetic field 
along the z axis (everywhere below we take H = 1 and 
the electron’s gyromagnetic ratio y e = 1, omitting the 
difference between the magnetic fields and the Larmor 
frequencies). Theoretical approach to the spin evolution 
in organic semiconductors customarily relies on the ap¬ 
proximation where the quantum hyperfine field given by 
the sum in Eq. 

N 

br = ^ ^ Qrfclrfc , (2) 

fc= 1 

is approximated as a static classical vector b r of ran¬ 
dom amplitude and direction, sampled from the Gaussian 
distribution with zero mean and the standard deviation 
equal to &hf = \\/Yhk a r fc> this approximation is justi¬ 
fied by the large number of nuclear spins coupled to the 
carrier spin at a given site (TV of order 10 or more)22r— . 
We also assume that the hyperfine fields at different sites 
are uncorrelated, so that hf = where 

a,/3 = x,y,z. 

When we consider a carrier with the initial spin along 
z, performing a random walk with trajectory r(t), the 
carrier’s spin fi(t) evolves according to the equation of 
motion 

R = b r(t) x fj,(t) = Cl(r(t))n(t) (3) 

where the matrix 

( o -K by \ 

fi(r)= h r 0 - fo r , (4) 

\-by b* 0 ) 



FIG. 2: (Color online) Diffusive random walk of a polaron 
(red circle) over d = 1 and d = 2 dimensional regular lattices 
of molecular sites (gray circles), (a) In the absence of exter¬ 
nal electric field, polaron hops between the nearest neighbor 
sites of d = 1 linear chain occur with the equal rates, zv/2, 
indicated by the blue arrows, (b) In the presence of external 
electric field along the linear chain, hopping rates Vf and Vb, 
respectively in the directions forward and backward to the 
electric field, are different; Vf > Vb. (c) With an external 
electric field along the x axis of a regular lattice in d = 2, 
hopping rates in the forward and backward directions to x, 
i/f and Vb, can depend on the electric field, while the rates in 
the perpendicular direction are the same as without the field, 
and equal to v/A each. 


describes the spin rotation taking place at the site r. 
Formal solution of this equation can be written in terms 
of the time-ordered exponent, 

fi(t) = T exp f df / fi(r(f , ))/x(0), 

J o 

with the initial condition /x(0) = z. The spin polarization 
is obtained by double averaging of the ^-component of 

rW, 


P(t ) = ((/^0)» = ((rexp/ dt'n(r(t'))\ ) 

W . n / hf / r 


(5) 


where (-) r w denotes averaging over the random walk tra¬ 
jectories, (-)hf denotes averaging over the local hyperfine 
fields, and the indices zz denote that we need to take the 
zz entry of the matrix which results after the averaging 
of the time-ordered matrix exponent. 

Without returns, all spin rotations at different sites 
would be independent and uncorrelated, leading to ex¬ 
ponential decay of the polarization as a function of time 
(the motional narrowing regime). In the presence of the 
returns, the rotations at different moments of time are 
correlated, and the polarization decay accelerates. In 
d = 1, the number of returns of the charge carrier to a 
given site after n hops is 0(n 1//2 ), whereas in d = 2 and 
d = 3 this number is 0(\n n) and 0( l), respectively^ 
Therefore one should expect that the influence of the re¬ 
turns is strong in d = 1, modest in d = 2, and weak in 
d = 3 dimensions; below we concentrate mainly on the 
d = 1 case where the effect is strongest. 
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The exact solution of Eq. (0 in the presence of returns 
is not available. In order to approach the problem, ev¬ 
erywhere below we employ the fact that the hopping is 
much faster than rotation in the hyperfine field, so the 
parameter 77 = bhf/v is small. Indeed, the value of frhf is 
typically of order of 100 MHz, while the average carrier 
hopping rate v is about 1-100 GHz, so 


77 = frhf/z' ~ 0.1 — 0.001 -C 1. 


We can calculate P(t) via the cumulant expansion in 
terms of the small parameter 77 : 


P(t) = 


Tex p/ dt'Cl(r(t')) 

Jo 


= exp I ^2 Knit )), 


where K n (t) is proportional to 77 ™. The odd cumulants 
vanish (since the local Gaussian distributions of b r have 
zero mean), the first non-vanishing cumulant K 2 (t) will 
determine the behavior of P(t), at least at short times. 
We will demonstrate the accuracy of this approach by 
comparing the analytically calculated K 2 (t) with the re¬ 
sults of the direct numerical simulations. 

The numerical simulations are even more important 
for studying the time-integrated polarization er(r). We 
are not aware of any analytical theory for this quantity 
which would provide insights and guide our investiga¬ 
tion. Thus, we rely solely on the numerical results, which 
demonstrate surprising and interesting universal features 
of er(r). 

One could do numerical simulations by Monte-Carlo 
sampling of the random-walk trajectories and the distri¬ 
butions of the local fields, thus calculating the average 
in Eq. 0 directly. However, our results show that the 
statistical error is quite large, so instead we employ the 
approach based on the Liouville equation. 

We describe the carrier spin via its density matrix 
Pr(t) = \ {q r (t) +m r (t)<r), where q T (t) is the probability 
to find the carrier at site r at time t, and m r (t) is its 
spin polarization; er is the vector of Pauli matrices. The 
carrier dynamics obeys the master equation 

= T l Wr '’ Tqr> (*) ~ W T,r'q r (t )], ( 6 ) 

r' 

where W r r > is the hopping rate from site r to r'. At the 
same time, the spin polarization follows the generalized 
drift-diffusion equation, 


dm r 

dt 


£[W r , r mi,(t) — W r>r 'in r (t)] + b r x m r . (7) 


For the charge carrier initially injected at the site r = 0 
in the spin-up state, the initial conditions correspond to 
q r ( 0 ) = S T fl and m r ( 0 ) = <5 r ,om(0) with m( 0 ) = ( 0 , 0 , 1 ) 
(directed along the z-axis). The solution m r (t) of Eq. 0 
includes averaging over the random-walk trajectories of 
the duration t, but the hyperfine fields at each site are 
taken as having some specific directions and amplitudes, 


i.e. the set {b r } of the local fields is fixed. Averaging 
over the local hyperfine fields is performed via Monte- 
Carlo sampling of {£» r } at each site from the Gaussian 
distribution A/"(b) = (27r6hf) -3 ^ 2 exp(— |b| 2 /26j( f ). 

In particular, in this way we determine the time- 
integrated spin polarization at a given location, 

/>00 

a(r) = v / dt{m z r (t)) h f, 

Jo 

which plays an important role in the spin transport mea¬ 
surements. 

The effect of the external magnetic field is included 
into our model by simply adding the external field B to 
the local hyperfine fields. The external electric field E 
applied along the x-axis (Figs. [T] and 0 is taken into ac¬ 
count by modifying the hopping rates: the hops along 
the field are more probable than backwards. Below, we 
assume, in the spirit of the Miller-Abrahams theory^ 
that the backward hopping rate Vf, (upwards in the elec¬ 
tric field potential) is exponentially suppressed in com¬ 
parison with the forward hopping rate Vf (downwards 
in the electric potential), i.e. Vb = {v/2) exp(— e), where 
e = eEa/kgT with the Boltzmann constant kg and tem¬ 
perature T, and eEa is the electric potential difference 
between two neighboring sites, while the forward-hopping 
rate remains unchanged, Vf = v/2. In particular, for 
d = 1 , we have 

W r y = {v/2)6 r y-i + {v/2)5 r y +1 for E = 0, ( 8 ) 

Wry = {v/2) S r y- 1 + (i'e -e /2) S r y +1 for E^ 0 . 

III. SPIN RELAXATION IN A LATERAL SPIN 
VALVE 

We neglect the effect of injector/detector electrodes, 
assuming insignificant tunneling between the leads and 
the semiconductor. For the lateral spin valve this im¬ 
plies unbounded diffusion over an infinite chain. For nu¬ 
merical simulations, we used a long chain with periodic 
or reflecting boundary conditions; the length was large 
enough to ensure vanishing population near the ends at 
all times. We also excluded from consideration the ad¬ 
ditional spin relaxation which is possible at the interface 
between a ferromagnetic electrode and an organic active 
layer— 


A. P(t) and a(r) in the absence of external fields 

Important insights about the short-time behavior of 
P{t) can be obtained analytically, using the lowest or¬ 
ders of the cumulant expansion in terms of the small 
parameter 77 = b^f/v, 

P{t) = exp ( ^2 K n{t) ) ~ exp [I< 2 
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renormalized time, z 


FIG. 3: (Color online) Logarithmic plot of P(t) in d = 1, 
for values of g = fchf/t' ranging from 0.002 to 0.05. When 
plotted against the renormalized time, r = (nt)g 4 ^ 3 , all data 
points fall on a single scaling curve. The resulting scaling 
curve is well fitted by exp(—8 t 3 ' 2 /3\/27r) for small r and 
0.827 exp(—1.33T 3 / 4 ) for larger t. Normal plot of the scal¬ 
ing curve, F(t), is shown in the inset. 


For a random walk over an infinite chain the first two 
non-vanishing cumulants, K% and A' 4 , are calculated in 
Appendix [B] The second cumulant K 2 is determined by 
the two-time correlation function, ((D(r(fi)) D(r(t 2 )))). 
From the large-t asymptotics of this correlation function 
we find: 




41 lL , 

\f 2 wv J J \Jt\ —12 3\/27r 

0 0 


(9) 

The next non-vanishing cumulant, A' 4 , is determined by 
the 4-th order correlation function of the process D(r(f)). 
Our calculations in Appendix [B] show that it has small 
numerical prefactor, K±{t) cs 0.01 • ? 7 4 (z 4 ) 3 , so that this 
cumulant becomes comparable to I \2 only at rather long 
times. Thus, K-iit) dominates the polarization decay at 
small times: 


2 

P(() » e*’«> = e-W')", = ( 10 ) 

Excellent accuracy of this scaling at short times is seen 
from comparison with the direct numerical simulations 
in Fig. [3] The similar decay law, P(t) ~ exp(—t 3 / 2 ), 
has been obtained in earlier studies, which assumed the 
single-axis local hyperfine fields (directed along the 2 - 
axis) 22 r— , or the hyperfine fields with fixed amplitude 
randomly distributed in the x-y planed. Our results 
confirm this decay law for the hyperfine fields distributed 
isotropically in space, and show that this feature holds for 
a very wide range of problems related to the spin decay 


during ID diffusion; we will also see the same decay law 
below, for the vertical spin valve case. 

More importantly, we notice that A 2 (f) and K 4 .it) de¬ 
pend only on the single dimensionless renormalized time 
r = (ut)r] 4 ^ 3 . From Eqs. (Q and ([5]) one can see that all 
cumulants, as well as P(t) itself, are the functions of two 
dimensionless quantities, (ut) and rj. However, our an¬ 
alytical and numerical studies evidence a much stronger 
result, that P(t) is a function of a single dimensionless 
quantity r. We performed a series of simulations for dif¬ 
ferent values of 77 , and Fig. [3] shows that all results fall 
on the same universal curve A(r), given in the inset of 
Fig. [3] This holds at all times we studied, even at large 
r, when the contribution from the high-order cumulants 
is important. 

Even more, we see that the same scaling holds when we 
consider the polarization decay at finite magnetic fields, 
finite electric fields, as well as for the case of the vertical 
spin valve (both without fields and with external mag¬ 
netic and/or electric fields). Thus, it is highly likely that 
the renormalized time r represents a universal feature of 
the spin decay for d = 1 random walk. Understanding of 
this remarkable scaling, as far as we know, is lacking. 

Another interesting feature, seen from Fig. 0 is the 
decay of P(t) at long times, which has a stretched- 
exponential form Pit) ~ exp(— ag{vt) 3 ^ 4 ), with a = 
1.33; we checked that this form remains very accurate all 
the way down to P(t) ~ 10 -12 . Again, to our knowledge, 
the reasons for this behavior are not understood yet. 

Equally interesting is the behavior of the time- 
integrated polarization cr(r). We are not aware of any 
analytical theory, which would provide insights in the be¬ 
havior of this quantity and guide our simulations. Thus, 
we rely solely on the numerical results. 

Our numerical simulations show that in the whole 
range of parameters a (r) has the exponential form, 
a(r) = lsexp(-\r\/l s ). Without external fields, the 
time-integrated polarization precisely follows the scaling 
law cr(r) = r]~ 2 ^ 3 G(rri 2 ^ 3 ), where the scaling function 
G(w) = 0.68 exp (—1.47|tu|) is obtained from numerical 
fitting. The origin of this scaling, as well as the origin of 
the exponential dependence of <r(r), are not clear. 

Note that the exponential decay of a(r) in the case 
of ID persistent diffusion is not trivial. If multiple re¬ 
turns were negligible (as in transient diffusion in 3D), 
P{t) would decay exponentially with the decay time 
ts, tr; the space- and time-integrated polarizations then 
would be related by the simple convolution^, cr(r) = 
v / Q °° dtP(t)q r (t), where q r (t) is the probability to find 
the carrier at the site r at time t, see Eq. ©• This would 
lead to the exponential decay a(r) = exp (— |r"|/Zsr,tr), 
with the well-known diffusion relation ls,tr = y/Dil, tr, 
where D is the diffusion coefficient (for a random walk 
on a d-dimensional lattice, D = va 2 /2d. where a is the 
average distance between the sites). 

However, in our case, where the returns are crucial, 
er(r) and P(t) are not related in such a simple way, and it 
is not even clear whether such a relation exists. Indeed, if 
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FIG. 4: (Color online) Magnetic field dependence of spin re¬ 
laxation for a diffusion in d = 1 with r/ = 0.01. Red, green, 
and blue plots correspond to the magnetic field values B = 0, 
0.5 bhf, and 2bhf, respectively, (a) P(t) is plotted against 
r = the red curve is the same as the one in the 

inset of Fig. [3] Inset: the same plot in the log-scale, (b) cr(r) 
is plotted versus r. Black curves are our exponential fits, (c) 
The normalized spin diffusion length, ls(B)/ls( 0), is plotted 
vs B/bhf- Dashed line indicates ls(B)/ls( 0) without returns, 
plotted from Eq. (1 A 131) . 


cumulants K 2 for zero magnetic field and A',f for finite 
magnetic field. This term induces a cutoff for t > B~ 1 , 
reducing the integral significantly. This is somewhat sim¬ 
ilar to motional narrowing: because of the external mag¬ 
netic field, the transversal components of the total field 
seen by /i, average out on the timescale B _1 . 

Comparing Eqs. ga and m, one can see that the 
cutoff induced by the external magnetic field becomes 
important for B ~ tg 1 ~ r / 1 which is even smaller 
than &hf- Our numerical results (Fig. [Tj clearly verify 
this sensitivity already at very low fields. This behavior 
is in striking contrast with the transient diffusion in 3D, 
where the spin relaxation time would scale as ts,t r oc 
(l + (B/i/) 2 ), meaning that the magnetic field effects 
would be visible only at very large fields B ~ v &hf- 

The time-integrated spin polarization cr(r) also ex¬ 
hibits strong sensitivity to the external magnetic field. It 
still has exponential form, cr(r) oc exp [— \r\/ls(B)] , but 
the spin decay length sensitively depends on B ; Fig.QJb) 
illustrates this dependence for 77 = 0.01. The magnetic 
field dependence of the (normalized) spin decay length, 
ls(B)/lg{ 0), is plotted in Fig. [2(c). Again, the analogy 
to the case of the transient 3D diffusion is superficial; 
this point is demonstrated in more detail in Appendix l~\1 
[see Eq. (IA13I) ]. where the transient diffusion case is an¬ 
alyzed, and its qualitative difference with our results for 
d = 1 are emphasized. 


C. Role of external electric field 


we used the same convolution of the non-exponential P(t) 
with q r (t) for ID diffusion, we would obtain clearly non¬ 
exponential decay law for cr(r). Thus, the origin of the 
exponential decay for d = 1 must be different from that 
of d = 3 case; this guess is supported by the qualitative 
difference in response of P(t) and cr(r') to the external 
magnetic and electric fields between the d = 1 and d = 3 
cases. 


B. Role of external magnetic field 

It has been noticed previously^ that in low dimen¬ 
sional diffusion, in the presence of multiple returns, P(f) 
is very sensitive to the external magnetic field. For the 
magnetic field B, the cumulant expansion of P(t ) can 
be carried out after applying the rotating-frame trans¬ 
formation, fi(t) —> exp(f flB)/x(t), where Qb is the skew- 
symmetric matrix formed of B (see Appendix iBl). Taking 
the external field as directed along the z-axis, we find the 
second cumulant 


t t\ 

Kz(t) = -2b'l f J dti J dt 2 

0 0 


cos(B[<i - t 2 \) 
sj2Tn;(ti - t 2 ) 


( 11 ) 


It is instructive to compare Eqs. m and ([9l): the cosine 
term in the integrand is the only difference between the 


If a drive voltage is applied to the spin valve (Fig. [T|), 
the resulting electric field E = Ex leads to a change 
in the hopping rates along and against the field di¬ 
rection (forward and backward hopping rates Vf and 
zz{,, see Fig. 0(b) ). Utilizing the Miller-Abrahams hop¬ 
ping models we take Vb = iy/2) exp(—e), where e = 
eEa/ksT , and Vf = v/2 (independent of E). Overall, 
this would lead to slower motion of the carrier, imply¬ 
ing slower changes of the random hyperfine field acting 
on it, and therefore (as it happens in the motional nar¬ 
rowing scenario) would produce faster decay of P(t) with 
increasing e. In the regime of 3D transient diffusion (see 
Appendix O for details) this is the most important ef¬ 
fect: P(t) would decay exponentially, and the decay time 
would decrease as (uf + Vb)/v- 

However, in the persistent-diffusion regime, the returns 
are important; besides inducing faster hopping, the elec¬ 
tric field also changes the statistics of the returns. By 
making the forward hops more probable, the probability 
of the returns is decreased, thus profoundly affecting the 
carrier spin relaxation. As above, we calculate the short- 
time behavior of P(t) using the cumulant expansion, and 
the second cumulant in the presence of finite electric field 
is 


Ff (f) = -6 2 f J dtiJdt 2 


e _ (v /F /^- \AU 2 [*1 -*2] 
{yfVbY^y/nit 1 - 1 2 ) 


( 12 ) 


0 


0 
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changes with electric field. 


IV. SPIN RELAXATION IN A VERTICAL SPIN 
VALVE 

The geometry of the vertical spin valve, Fig. [Q)b), sug¬ 
gests diffusion over a linear chain of a finite length L. 
Neglecting the back-tunneling into the electrodes, as it 
often happens in experiments, we obtain the perfectly re¬ 
flecting boundaries. The spin relaxation now depends on 
the length of the system: P(t) = P(t,L), a(r) = a(r,L). 
Another feature of this geometry is that, instead of the 
whole function a(r,L), one is interested in its value at 
the detection electrode, cr*(L) = <r(L,L). 


A. P(t) and in the absence of external fields 


FIG. 5: (Color online) Electric field dependence of spin relax¬ 
ation in d = 1 with r/ = 0.01, at zero magnetic field. With 
typical a as 1 nm and T as 100 K, parameter E = eEa/ksT 
corresponds to the electric field E as 8.6 e mV/nm. The red, 
dark green, and violet plots correspond to e = 0, 0.15, and 
0.3, respectively, (a) P(t) is plotted versus r = (vt) jj 4 / 3 ; the 
red curve is the same as in Fig. [H)a). Inset: the same in the 
log-scale, (b) a(r) is plotted versus r. Black lines are our ex¬ 
ponential fits, (c) Normalized spin diffusion length ls{e)/ls( 0) 
versus e (orange). The dashed line indicates the same depen¬ 
dence without multiple returns, plotted from Eq. (1A16I) . 

Here the exponent shows that the electric field prevents 
multiple returns. Because of the exponent, the integral 
in Eq. m mainly decreases with increasing E (see Ap¬ 
pendix [B] for more details). This means that the spin 
relaxation slows down with increasing electric field, at 
least at relatively short times. The expected effect is 
seen in Fig. [5ja) , which demonstrates simulation results 
for P(t) for three different values of E. 

Now we turn to the time-integrated spin polarization, 
<r(r). Our numerical analysis shows that, like in all cases 
above, in the presence of finite electric field <r(r) has ex¬ 
ponential form, a(r) oc exp [— r/ls(E)\; see an example 
in Fig. [U)b) for r/ = 0.01. The curves are not symmet¬ 
ric with respect to r = 0 , which is an obvious result of 
the drift induced by the electric field. Also note that, in 
contrast with two previous cases (no external fields and 
external magnetic field), the value of <r(0), which sets 
the magnitude scale for the whole curve, is not propor¬ 
tional to ls{E) anymore, and decreases with increasing 
E, although ls(E) itself increases. 

The spin decay length is very sensitive to the electric 
field: in Fig. [Sjc) we show the dependence of ls(E) on 
the normalized electric field e = eEa/ksT for rj = 0.01, 
where Zg(0) ~ 14.9. This is to be contrasted with the sim¬ 
ilar dependence for a transiently diffusing carrier in d = 3 
with the same ls{0 ), which is shown in the same graph for 
comparison: without returns, ls{E) shows much weaker 


Let a carrier be implanted at the boundary site, r = 1, 
of a linear chain of finite length L. After n hops it will 
diffusively cover the distance ~ y/n. With the hopping 
rate v one has n ss vt, so that for relatively short times, 
vt < L 2 , the boundary at r = L will not affect the spin re¬ 
laxation noticeably. Therefore, for relatively short times, 
P(t) can be found by considering a carrier diffusing over 
the semi-infinite chain, r = 1 , 2 ,..., with the reflecting 
boundary at r = 1. Calculation carried out in Appendix 
[B]for this case gives the second cumulant function, 

K>{t)^-^={vtf/\ (13) 

which differs from K 2 of the infinite chain, Eq. m, only 
by the factor y/2. The ensuing short-time superexponen¬ 
tial dependence of P(t) is confirmed in our simulations, 
see Fig. [ 6 ] 

At longer times, the influence of both boundaries be¬ 
comes noticeable, and the decay of the spin polarization 
will depend on L. This dependence can be guessed using 
the results for unbounded diffusion given above. 

As we have seen in the previous Section, the character¬ 
istic length of the spin relaxation is Is, and scales with 
= bhf/v as t; - 2 / 3 . Based on this fact and on Eq. (flfil) . 
we can guess that the spin polarization should depend on 
the dimensionless time r = (z/t )^ 4 / 3 and on the dimen¬ 
sionless length A = Li r / 2 ^ 3 . Our numerical simulations 
confirm the expected scaling law P(t,L ) = F(t, A), and 
provide the most notable features of the scaling function 
F(t, A). In Fig. [G] (a)-(c), we demonstrate the scaling 
by plotting P(t,L) as a function of r and A, for three 
different values of A and twelve different values of ij. 

These figures also show that the scaling function has 
the superexponential form F(t, A) = exp(— 8t 3 / 2 /3^tt) 
at small times (r < 0.5), in accordance with Eq. (fl3l) . and 
is independent of the normalized chain length A. At large 
times, the scaling function is accurately described by the 
stretched exponential, F(t, A) ~ exp(— a(A)r^ A )), with 
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FIG. 6: (Color online) Spin relaxation in a finite chain of 
length L. When plotted versus r = ( vt ) r; 4 / 3 and A = Lr/ 2 ^ 3 , 
including (a) A = 0.75, (b) A = 1.25, and (c) A = 2, all P(t, L)- 
points fall on universal curves. Cyan lines are stretched ex¬ 
ponential fits, P(t,L) ~ exp (—at 13 ), with (a) /3 = 0.59, (b) 
/I = 0.61, and (c) /3 = 0.64. Magenta lines are our theoret¬ 
ical fits with exp(—8r 3 ^ 2 /3y / 7r). (d) The resulting universal 
scaling function, F(t, A), is plotted for the above values of A. 
The dotted line illustrates saturation of F(r , A) for A —» oo. 


FIG. 7: (Color online) Spin polarization at the detection elec¬ 
trode, a*(L), for a system with r/ = 0.01. (a) Log-plot of 

a *(-L) in zero electric field and different magnetic fields: B = 0 
(black), 0.56hf (red), 2fehf (green), (c) Log-plot of <t*(L) in 
three different electric fields, e = 0 (black; the same as in 
(a)), e = 0.03 (blue), e = 0.06 (magenta), in zero magnetic 
field. The black and cyan lines in (a) and (c) are our ex¬ 
ponential fits, (b) and (d): Dependence of the (normalized) 
diffusion length, ls/ls{0 ), on external magnetic and electric 
fields, respectively. 


the parameters which depend on the chain length; specif¬ 
ically, /3(0.75) = 0.59, /3(1.25) = 0.61, and /3(2) = 0.64 in 
Figs. [ 6 ] (a)-(c), respectively. As the length of the chain 
increases, the exponent /3(A) increases, saturating at the 
value j3 = 0.75 for very large A that corresponds to the 
unbounded diffusion [see Fig. [U(d)]. 

In a similar way, we numerically verify the existence of 
scaling for the time-integrated spin polarization. As we 
expect, all length scales are scaled by the factor r/~ 2 / 3 , 
so that a(r,L) = r]~ 2 ^ 3 G[rr] 2 ^ 3 , Lr] 2 / 3 ), in analogy with 
the case of unbounded diffusion in the lateral spin valve. 
Correspondingly, for the spin polarization <r*(L), ob¬ 
served at the detector electrode, we have the scaling 
<t*(L) = ?? _ 2 / 3 G*(L?y 2 / 3 ). The numerical fitting shows 
that the scaling function is very accurately described as 
exponential, G*(w) = 2.712exp(—1.475w). This corre¬ 
sponds to the dependence, cr*(L) = 4/5 exp (—L/ls), with 
the diffusion length, lg = 0.678 r]~ 2 ^ 3 , which is nearly 
identical to the one found for the unbounded diffusion in 
the lateral spin valve. 


B. Role of external magnetic and electric fields 

In analogy with the case of the lateral spin valve, we 
studied the influence of the magnetic field Bz, along the 
z-axis, and of the electric field E~k directed along the 
£-axis, on the spin relaxation. 

We are primarily interested in the behavior of er*(L). 
Our simulations (Fig. Q show that it remains essentially 
exponential with L , in both external magnetic and elec¬ 


tric fields. Thus, the field dependence of cr*(L) is fully 
encompassed by the field dependence of the spin decay 
length on B and e = eEa/ksT, shown in Figs. Ob) and 
(d). As before, we see that the returns lead to strong sen¬ 
sitivity of Is to the external fields, and this dependence 
is very close to its analog established in the previous Sec¬ 
tion. The curves in Figs. Ob) and (d) appear to resemble 
the ones in Fig. O c ) and Fig. O c )- 

V. SPIN RELAXATION IN d = 2 

Above we focused on the analysis of the spin relaxation 
of a diffusing carrier for d = 1 , mainly because the effect 
of multiple returns is strongest in this case. Meanwhile, 
it is rather straightforward to extend our analysis to the 
d = 2 case. For a carrier performing a simple random 
walk over a regular lattice in d = 2 , the second cumu- 
lant function is calculated in Appendix [B] It suggests 
the short-time decay 

P(t) ~ exp[— 2rf‘(vt) , (14) 

where 7 = 5.243 is a numerical coefficient. We have 
checked numerically that this formula accurately de¬ 
scribes the spin relaxation down to rather small polar¬ 
ization values. Specifically, at P(t) ~ 0.05, the observed 
deviation from Eq. m was only about 0.002. Our sim¬ 
ulations have also shown that at longer times the de¬ 
cay slows down, closely resembling exponential. This is 
shown in Fig. [HKa), which illustrates the spin relaxation 
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(c), where we compare the results for the random walk 
in different dimensions for 77 = 0.025, where the spin 
diffusion length in the absence of the external fields is 
Is = 9.0. From the results of Section Hill one can see that 
the same zero-field spin diffusion length in d = 1 corre¬ 
sponds to 771 ss 0.021, whereas for d = 3, from the strong 
collision approximation (i.e., by neglecting multiple re¬ 
turns) we find the same zero-field spin diffusion length for 
773 « 0.032 [cf. Eq. (I A13 [) ] . Therefore, in Figs. 0b) and 
(c) we compare the field dependence of the normalized 
spin diffusion length, ls{B)/ls and ls(E)/ls, for d = 1 
with 771 = 0.021, d = 2 with 77 = 0.025, and d = 3 
with 773 = 0.032. Evidently, more frequently occurring 
multiple returns lead to stronger growth of spin diffusion 
length with external fields. 


FIG. 8 : (Color online) Spin relaxation of a carrier diffusing 
over d = 2 regular lattice, with 77 = bhf/t' = 0.025. (a) Log- 
plot of P(t) against vt (black). The cyan dashed curve is 
plotted from Eq. m- The magenta dashed line is our expo¬ 
nential fit. (b) Magnetic field dependence of the (normalized) 
spin diffusion length; ls/ls{ 0) is plotted vs B for a diffusion 
in d — 1 with 771 = 0.021 (orange), d = 2 with 77 = 0.025 
(black), and strong-collision approximation to d = 3 with 
773 = 0.032 (dashed line), (c) Electric field dependence of the 
spin diffusion length. ls/ls( 0) is plotted vs e = eEa/ksT, 
for a diffusion in d = 1 with 771 = 0.021 (orange), d = 2 
with 77 = 0.025 (black), and d = 3 with 773 = 0.032 in strong- 
collision approximation (dashed line). 

in a system with 77 = 0.025. Thus, in d = 2 the effect 
of multiple returns is quite noticeable. Meanwhile, it is 
rather easy to check that in d = 3 the polarization decay 
is exponential at virtually all times. 

To understand the spin-transport relaxation in d = 2, 
we consider a lateral spin valve device similar to that 
of Fig. 0a), where the carriers hop between the sites 
of d = 2 regular lattice located in the x-z plane [see 
Fig. He)], so that each site is characterized by the radius 
vector r = (x,z), with x.z = 0,±1,±2,... Assuming 
that the size of the organic layer in z-direction is much 
larger than the spin diffusion length, we characterize the 
spin-transport relaxation by the quantity, 

poo 

d{x)= ^2 <j( d c ) = v X! / dt(m^.(t)) h {, (15) 

r =(*.*) r = (x,z) 0 

x=L x=L 

which is the total time-integrated spin polarization that 
has reached the detection electrode at x = L. Here m T (t) 
is the solution of drift-diffusion equation 0 for the 2D 
regular lattice; we solve this equation numerically, with 
the initial condition m r (0) = <5 r) om(0). The influence 
of the external magnetic and electric fields is taken into 
account as described above. 

Our results show that the decay of U 2 d{x ) is exponen¬ 
tial, both with and without the external fields. The spin 
diffusion length lg is also rather sensitive to the exter¬ 
nal fields, as a result of the returns in the course of the 
2D random walk. This is illustrated in Figs. 0b) and 


VI. CONCLUSION 

We have investigated spin relaxation of a carrier per¬ 
forming a random walk on a lattice, with random mag¬ 
netic fields at each site. This models spin relaxation 
in organic semiconductors, where the charge transport 
between the 7r-conjugated segments of molecules is in¬ 
coherent, and where the carrier spin interacts with the 
hydrogen nuclear spins surrounding a segment. Due to 
relatively large number of surrounding nuclear spins and 
their slow dynamics, the on-site random magnetic fields 
are taken as static, and sampled from the Gaussian dis¬ 
tribution. The width of the distribution reflects the 
strength of hyperfine coupling at each site. For the in¬ 
coherent hopping we assumed a random walk with the 
constant transition rates. 

To understand the effect of multiple self-intersections 
of random walks, we have focused on the motion in 
d = 1 dimensional infinite linear chain. A superexpo¬ 
nential short-time decay of spin polarization, P(t) ~ 
exp(-aif 3/2 ), was found analytically from the curnu- 
lant expansion. The numerical simulations confirmed 
the superexponential dependence, and showed that it 
changes to a stretched exponential at longer times. We 
have also analyzed the spin relaxation of a carrier dif¬ 
fusing over a linear chain of finite length, and estab¬ 
lished that the short-time relaxation is somewhat faster, 
P(t) ~ exp(—\/2 q!i f 3 / 2 ), whereas the stretched exponen¬ 
tial long-time decay becomes slower. As a consequence 
of the multiple returns, in all these cases P(t) is highly 
sensitive to the external magnetic and electric fields. For 
a diffusion over d = 2 regular lattice, we have found 
the short-time behavior P(t) ~ exp(—Inf), smoothly 
crossing over to the exponential long-time decay. 

Due to its relevance for the spin transport experiments, 
we have investigated the spin relaxation in the space do¬ 
main, considering the time-integrated spin polarization 
cr(r). It was demonstrated that, despite the strongly non¬ 
exponential decay in time, the spin relaxation in space is 
essentially exponential, superficially similar to that of a 
carrier diffusing in d = 3. However, diffusion in lower di- 
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mensions shows much stronger sensitivity to the external 
electric and magnetic fields. 
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Appendix A 

In this Appendix we calculate spin relaxation for a 
carrier diffusing transiently, i.e., when self-intersections 
of its random-walk trajectories are negligible, and the 
strong-collision approximation^! is valid. Our starting 
point is Eq. © of the main text. Let r (t) be a random 
walk trajectory that starts at r 0 = 0 and passes through 
n sites, iq, ..,r„. Then the time-ordered exponent can be 
written as 


Next we want to average Eq. (IA2I) over the waiting 
time distributions and random-walk trajectories. Be¬ 
cause of the absence of returns, the latter reduces to a 
summation over all n, whereas the former can be done by 
integrating f(r) with the waiting-time distribution func¬ 
tion, v e~ VT . Hence from Eqs. © and m we get: 


p (t) = E I d,To ■ ■ ■ j dr n Y[f {rj)ve 

3=0 


n—0 


0 0 

n —1 \ / n 

t - e T k ) -0 ( t - E Tk 

k—0 ' ' k=0 


(A5) 


Here, the difference of 0-functions guarantees that at 
time t the walker has performed exactly n steps, so that 
J2kZa T k < t < J2k=o Tk ~ Using the integral represen¬ 
tation, 0(x) = f[e lzx /(z — ie)]dz/(2ni), we reduce Eq. 
(IA5I) to 


P(t) = 


I" dz e izvt u{\)-u{l+iz) 
J 27 n z — ie 1 — u{ 1 + iz) ’ 


— OO 


(A 6 ) 


where u(y) = v J^°drf(T)e yl/T . Going back to the defi¬ 
nition of /(t) ant taking the r-integral we get: 


Texp f di / n(r(t')) = e T “ n(p "> • • • e T ° n(ro) , (Al) 
Jo 

where r k are the waiting time at r*,, k = 0,1,.., n. In 
the absence of self-intersections all rq are different, and 
the average of Eq. EU) over the local hyperfine frequen¬ 
cies is a product of exponents, averaged over the inde¬ 
pendent Gaussian distributions of hyperfine frequencies 
{b rfc }, namely 


u(y) 


I [fir 6 — (ot) 2 + (Cz+/3) 2 

vJ ^ (2n)i {yriY + tf + Q + ^ + Pr 


(A7) 

where /? = B/b hf. Exact evaluation of this integral yields 
u(y) in terms of the error function. However, we are 
interested in small values of r] = bhf/v, where the most 
relevant pole of Eq. (lA6l) . given by u( 1 + iz 0 ) = 1 , is 
located at \zo\ <C 1. Therefore zq can be found from the 
small- 77 (77 -C /3 _1 , 1) expansion of Eq. (I A 711 for \y\ ~ 1, 


T exp / dt'fl(r{t’)) 


hf 


=n 


(A2) 


k-0 


where /(r) is defined by (e r ° b ){b} = i/(r), with fib 
being the skew-symmetric matrix formed of b. By aver¬ 
aging 


U{y) ^y { 1 (y/ 77)2 + / 3 2 + 1 )> (A8) 

which in fact provides a good approximation for any /3. 
Form Eq. (IA 8 I) we find zq ~ 2i/{r]~ 2 + (3 2 + 1), yielding 

P(t) * e-^'s , 4 = (T 7 " 2 + /3 2 + l)/2v. (A9) 


e rnb = i +sin(|b|r)^ + 2 sin 2 (|b|r/ 2 )^, (A3) 

over the Gaussian distribution of b with zero mean and 
standard deviation, b^ f, one gets: 

f( T ) = ^ + ^(! - & hU 2 ) exp(- 6 hU 2 / 2 )- ( A4 ) 


The spacial dependence of spin polarization is given by 

°? 00 °? °? 

<r{ r) = v dt y^Q w (r) dr 0 ■ ■ ■/dr n ]~[ f{Tj)ve~ VTi 


3 =0 


0 ( t - r fc j -0 f t - r k 

k—0 ' ' k—0 


, (A10) 


When a longitudinal magnetic field B = Bz is applied, 
Eq. (1A3I) should be averaged over a Gaussian distribution 
of b with the mean, (b) = B z. Even though this leads 
to a non-diagonal matrix (e r ° b ){b}j it remains block- 
diagonal, so that Eq. (IA2I) holds with a modified /(t). 
Note in passing that /(r) is a typical example of the 
static Kubo-Toyabe relaxation function ! 41 > 42 


where Q n (r) is the probability that the random walker is 
at r after n steps. Calculating the integrals in Eq. (IA10I) 
is easy by taking first the f-integral. Further, introducing 
u = v 2 j^ 1 dTTf(r)e~ l ' T and using Eq. (jA 8 [) . we find: 

OO OO 

cr(r) = «E <2™(r)[u(l)] n ^ ^2Qn{r)e~ n/vt ' s . (All) 

n =0 n —0 
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The probability Q n { r) is related to the solution of Eq. © 
as q r (t) = e~ ut Q n (r)(ut) n /n\. This can be used 

in Eq. (1A11I) to express <r(r) in terms of q r (t ): 

/*00 

er(r) ~ v / dtq r (t)P(t). (A12) 

Jo 

The large -r behavior of <r(r) follows from that of 
q r (t). Namely, for a simple random walk on 
a d-dimensional regular lattice one has g r (£) = 

(27 Tvt/d)~ d / 2 exp[—d|r| 2 /( 2 ^t)], leading to the exponen¬ 
tial decay, er(r) oc exp(— |r|/Z^), with 

= + (A13) 

In an external electric field, E = Ex, hopping rates 
along x are changed. This leads to the drift along x, 
and also modifies the diffusion and the waiting-time dis¬ 
tribution. While P(t) is affected only because of the 
change in waiting-time distribution, q r (t ) and conse¬ 
quently <r(r) are sensitive to the drift and diffusion. As¬ 
suming a random walk over a d-dimensional regular lat¬ 
tice, and for the hopping model considered in the main 
text, the hopping rates forward and backward to x are 
Vf = v/2d and v b = ve~ e /2d 1 where e = eEa/ksT , 


whereas in perpendicular directions hopping rates are 
v/2d. From the corresponding waiting-time distribution 
function, ve~ VT with v = Vf + v b + ^(d — l)/d, one finds 
P(t) ~ exp[—£/£^(e)] with the electric-field dependent 
spin relaxation time, 

, 1 2 d—l + e~ £ p 2 + l d 

s 2vq 2 2d v 2d—l + e~ £ 

(AH) 

The drift-diffusion equation © in d = 3 dimensions has 
the solution, 

q r {t)= j e~ vt I x ( 2 y/vf v b t ) I y ( vt/d ) I z (vt/d ), 

(A15) 

where r = ( x,y,z ) with x,y,z = 0 ,± 1 ,± 2 ,.., and I a is 
the modified Bessel function of order a. For a spin valve 
similar to those illustrated in Fig. ©we evaluate the 
quantity, a(x) = „ <r(r) (in lower dimensions, one or 

both of last terms in Eq. (IA15I) should be eliminated, and 
the sum for cr(x) should be changed correspondingly). 
After taking this sum, the integral Eq. (IA12I) reduced 
to the Laplace transform for a modified Bessel function, 
yielding a(x) oc exp[— x/ls{e)\ for all x > 0 , where 


l s(e) = 


In i(l + e- £ ) + 


PiCJs) 


+ 


(A16) 


i(l + e- E ) + 


d 

vt's^s) 


This dependence is plotted in Fig. M c) for a system with /3 = 0 and ls(0) = 14.9 (corresponding to q = 0.0336 in 
d = 1), and in Fig. ©c) for a system with /3 = 0 and l${ 0) = 9 (corresponding to q = 0.0322 in d = 3). 


Appendix B 

In this Appendix we calculate the d = 1 dimensional second cumulant functions, A' 2 , A'.f, Ajf, A'^, Eqs. m, 
mu), m , and m the fourth cumulant function K 4 , as well as the second cumulant for a simple random walk on a 

( 2 ) 

two-dimensional regular lattice, K\ . Basic ingredients of this calculation are the Markov property of random walk 
and its Greens function, G(r, r',t), which is the solution of corresponding random walk equation © with the initial 
condition, q r '(0) = S r y. The second cumulant function of Eq. © is defined by 

t ti 

K 2 (t ) =JdhJ d£ 2 «n(£ 1) n(£ 2 )))„, «n(tr) n(£ 2 )» = ((L2(r(£ 1 )) H(r(£ 2 ))) hf ) rw , (Bl) 

0 0 

where the average over random walk trajectories and locally Gaussian hyperfine frequencies is meant. Using the 
matrix form, Eq. ©, and the fact that the components of b r are delta-correlated, ( 6 “ = ^hA-.r'CG,/?, one 

easily takes the average over local frequencies, resulting in ((fi(£i) 0(£ 2 ))) = —1 ■ 2 (5 r (t 1 ), r (t a )) rw - F° r any type °f 

random walk, this can be expressed via the Greens function as follows: 

((0(£r) n(£ 2 ))) = ~i • 2 &h f G ( r ’ r > G “ ^)G( 0 , r, £ 2 ). 


(B2) 
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For the random walk on an infinite chain we have G(r,r',t) = e vt I\ r _ r i\(vt ), where I r (z) is the modified Bessel 
function of order i\ From Eq. (IB2I) we find, ((fl(fi) ft(t 2 ))) zz = — 2 b^ l{ e~ l '^ tlt2 ^o {y(t\ —f 2 )), which gives the second 
cumulant, 


Ut Z\ 

K 2 (t) = —2 if J dzi jdz 2 e~ ( - Zl ~ Z2 ' > Io(zi — z 2 ). (B3) 

o o 

In view of large 77, it is necessary to find the integral for large (ut) 1. Utilizing the large-z asymptote e~ z I r (z ) ~ 
(27T2) -1 / 2 exp(— r 2 /2z) in the integrand, we arrive at the result Eq. ([2]). 

As the odd cumulants are zero, the fourth cumulant function is expressed in terms of the four-time correlation 
function as follows: 


t t\ £2 £3 


+ K 4 (t) = jdt 1 Jdt 2 Jdt^jdU{{n{t l )Vt{t 2 )Vl(t 3 )Vl(t A ))) zz . 


(B4) 


0000 

From Eq. JI]) we find that the zz-component of the product fl(U) fl(t 2 ) f2(f 3 ) fl(t 4 ) is equal to 


h z h z 

b r(t 2 )°r(t 3 ) 






(B5) 


To find (fi(ti) fl(t 2 ) fl(t 3 ) tl(t 4 )) zz , we first average Eq. (1B5I) over the local hyperfine field distribution, then over the 
random walk trajectories. After the first averaging we get: 

(H(tl) Jl(t 2 ) ^(^4)}zz = &hf (4 ^r(ti),r(t 2 )^r(t 3 ),r(t 4 ) + 2 ^ r (t 1 ), r (t 3 ) < 5 r (t 2 ), r (t 4 ) + 4 (ir(ti),r(t 4 )^r(i 2 ),r(t 3 )) rw - (B6) 

This equation follows from the calculation of local field averages of the form, (b^b^b^b^) . For a ^ /3, this 
calculation is simple and gives b^ { <5 ri . r 2 <5r 3 , r4 - For a = /?, on the other hand, it results in the combination, 
^hf(^ri,r 2 ^r 3 ,r 4 + <5ri,r 3 (fr 2 ,r 4 + <5 ri ,r 4 d T2 ,r 3 )■ Note that contributions with r 3 = r 2 = r 3 = r 4 cancel out from this 
combination due to the Gaussian character of local frequency distributions. 

Next we average Eq. (IB 6 I) over the random walk trajectories. For a function of four coordinates, /, and times 
arranged as U > t 2 > t 3 > £ 4 , from the Markov property of random walk one generally has: 


(/(r(ti),r(t 2 ),r(t 3 ),r(t 4 ))) rw = ^ /(ri,r 2 ,r 3 ,r 4 )G(r 12 ,i 12 )G(r 23 ,t 23 )G(r 34 ,t 34 )G(r 4 ,t 4 ), (B7) 


where ry = — Yj and tij = U — tj. We apply Eq. (1B71) with the infinite-chain Greens function to each term of 

Eq. (fB 6 l) . and find the large- utij asymptotes of the resulting quantities: 

(B 8 ) 
(B9) 
(BIO) 


(^r(ti),r(t 2 )^r(t 3 ),r(t 4 )) rw G( 0 , tl 2 )G( 0 , £34) ~ 

(^r(ti),r(t 3 )^j-(t 2 ),r(t 4 )) rw = ^ ' G(t , t\ 2 )G(l *, t 23 )G(r, £ 34 ) ~ 


2lTU^/t 13 t 2A - t\ z 


(^r(ti),r(t 4 )^r(t 2 ),r(t 3 )) rw G(0, £23) ^ ( G(r, ti 2 )G(r : £ 34 ) — 


t 2 

c 23 


270/^44^23 - 

Equations (IB 6 I) and (IB8I) - (IB10I) define the integrand of Eq. (IB4I) . Using the asymptotic forms, we find: 

£ £1 £2 £3 

f dt 1 J dt 2 J dt 3 J di 4 ((fi(ti) fl(t 2 ) ft (£ 3 ) £1{U))) ZZ ~ ^if(utf. 


(Bll) 


0000 


This relation, together with Eq. (IB4D . leads to the result, 
Ki{t) ~ [( 57 t — 16)/97r]?7 4 (i/f) 3 . (B12) 


It is worthwhile to notice that Eqs. (IB9I) and (IB 101) 
make the non-Gaussian character of the stochastic pro¬ 
cess {!!(£)} explicit. Indeed, a Gaussian {fl(£)} would 
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entail a factorization of four-time functions, as it hap¬ 
pens in Eq. ( 11181 ) . whereas Eqs. m and mm do not 
satisfy this condition. 

Consider now the spin polarization decay in d = 1, in 
the presence of a magnetic field along z. This case can be 
described by adding in Eq. © the term B x m r , where 
B = B z. A straightforward evaluation of ((f2(fi) fl{t 2 ))) 
by repeating the steps that have led from Eq. m to 
Eq. (IB3I) is insufficient; magnetic field effects appear 
only in the fourth order, as a correction to K 4 of or¬ 
der ~ (Bbhi) ■ Rather, a systematic expansion of P(t) 
in powers of &hf can be achieved after performing the 
rotating-frame transformation, m r (f) = exp(t 0 ^) 11 ,.(t). 
Then n r (t) satisfies Eq. © with time-dependent local 
frequencies, b r (t) = exp (—t f7s)b r exp(f Clg), and ini¬ 
tial condition, n r (0) = d r om(0). Also, as exp(f does 
not change the 0 -component of the vector on which it 
acts, P(t) is expressed through n r (i) exactly in the way 
it was in terms of m,.(f). Averaging over the distribu¬ 
tion of local hyperfine fields now gives {{b(ti) b{t 2 ))) zz = 
-2 bl { cos(B[ii - ^])(<5r(i 1 ),r(t 2 )) rw , which leads to Eq. 
m after applying Eq. (IB2I) and large- {vt) expansion of 
the resulting Bessel function. 



FIG. 9: (Color online) exp[K 2 {t)\ is plotted against {vt) from 
Eq. (IB13I) with y = 0.05, for e = 0 (red), e = 0.5 (green), 
and e = 1 (blue). Open black circles are simulated points for 
P(t) with corresponding values of parameters. 

In the presence of electric field along x, the hopping 
rates forward and backward to x, Vf and 77 ,, are dif¬ 
ferent. In the paper we utilize vf = zx/2, and = 
{v/2)e~ E . The diffusion propagator is given by q r {t) = 
{vf/vb)i e~ ( ' Vf+Ub ' >t I r (2 y/VfVb t ) ■ This propagator leads 
to the second cumulant function, 


t 1 1 

Ki(t) = -2bl f JdtiJdt 2 e~ ( -‘' f+Vb)[tl ~ t ^Io (2y/vfv b [ti - t 2 }) ^ 
0 0 


bl c r r [* 1 —* 2 ] 

- - —jj / dti / dt 2 - =- 

( v fVb) 1/4 J J \J 7r(ti - t 2 ) 

00 v 


(B13) 


The ^-dependence of K 2 (t) is non-trivial because the prefactor in Eq. (IB 131) grows as exp(e/4), while the integral is 
suppressed with E. To clarify this dependence and to show that, for almost all t, the absolute value of the cumulant 
decreases with increasing E, we plot exp[/\|'(t)] with 77 = 0.05 in Fig. [9l for three different values of E. 

To find the cumulant function K% (t) for the case of a random walk over the semi-infinite chain, r = 1, 2,..., with the 
reflecting boundary at r = 1, we exploit its Greens function, G > (r,r',t) = e~ vt [/| r _ r /| {vt) + I r+r '-i{vt)\. Plugging 
this Greens function in Eq. (IB2I) and using the large-^f expansion of the modified Bessel function yields 


{(n(ti)n(( 2 ))) 



V y/ti — t 2 


Vh + 1 2 


(B14) 


By further integration we find the cumulant function Eq. m- 

The regular lattice in d = 2 can be described by the radius-vector, r = (x,y), with x,y = 0,±1,±2,..., 
and x, y along the lattice sides. The Greens function of the random walk over this lattice is G^ 2 ^(r,r',t) = 
(yt/2')I\y_ y i\ (vt/ 2). Further extension of Eqs. (IB1I) - (IB3I) to d = 2 is straightforward, leading to 


t 1 1 

K 2 \t) = -2 6 hf J dti J dt 2 - t 2 )/ 2 ) 

o o 


Ft 

—2 y 2 J dz{vt — z)c~ z Iq(z/2). 
o 


(B15) 


Combining numerical integration and asymptotic expansion of Iq(z) for large z, we find the large -0 expansion, 


/ > 1 1 

dz’{z — z’)e~ z Iq(z’ /2) = —zlw^z) — —\w{z) + </+ 0{1/ z), 7 « 5.243, ( ~ —0.264, (B16) 

' 7 r 27T 

0 


yielding 


K^it) ~ —2y 2 {vt) ln(7z/t)/-7r. 


(B17) 
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